Nonequilibrium noise in electrophoresis: the microion wind 



buropnya and Sriram RamaswamjQ 

Indian Institute of Science 
(Dated: July 2, 2009) 



ON 
O 
O 

(N 



(N 

S 1 

O 

£ 

i 

-t— > 



O 
o 



> 

ON 
NO 

q 
o 

On 
O 



X 



We analyze theoretically the dynamics of a single colloidal particle in an externally applied electric 
field. The thermal motions of microions lead to an anisotropic, nonequilibrium source of noise, pro- 
portional to the field, in the effective Langevin equation for the colloid. The fluctuation-dissipation 
ratio depends strongly on frequency, and the colloid if displaced from its steady-state position relaxes 
with a velocity not proportional to the gradient of the logarithm of the steady-state probability. 



PACS numbers: 82.70.Dd, 05.40.-a 



Colloids in electric fields present challenging problems 
in nonequilibrium statistical mechanics, as a result of the 
interplay of hydrodynamic and electrostatic interactions 
[Hj 0- Electrophoresis and aggregation into chains, crys- 
talline phases and fractals are just some of the diverse 
phenomena seen in these systems Q, However, a 

thorough theoretical understanding of these effects must 
begin with a single colloidal particle in an electric field. 
Important progress in this direction has been made by 
Squires 6] who has shown that particles stably suspended 
by the interplay of external fields, boundaries and hydro- 
dynamic flow at low Reynolds number respond to distur- 
bances as though they were in equilibrium in an effec- 
tive potential. The effective particle velocity divided by 
the Stokes drag is treated as a pseudoforce which sub- 
sumes all hydrodynamic interactions of the colloidal par- 
ticle and surrounding counterions with the wall. This 
pseudoforce, it turns out, can be written as the gradi- 
ent of a scalar 'pseudopotential'. A major success of this 
approach was the explanation of apparent like-charge at- 
traction in certain experiments j 01 ■ 

To extend this formalism beyond the average motion 
of the colloid, note that it is subject to two kinds of 
stochastic forces. One, which is taken into account in 
Q , is thermal noise corresponding to viscous dissipation, 
with variance proportional to temperature times Stokes 
drag. The other enters as follows: the local charge den- 
sity p fluctuates because of thermal motions of the mi- 
croions; in the presence of an imposed electric field E, 
this means the electric force density pE in the Stokes 
equation fluctuates as well. This results in a fluctuating 
contribution to the motion of the colloid, with variance 
proportional to temperature times the square of the elec- 
tric field and correlations controlled by microion motion 
Q. This nonequilibrium coloured noise, ignored so far as 
we know in earlier work [fj, is the main subject of this 
Letter. 

We summarize our main results before presenting de- 
tails of the work. We have studied two systems: (I) a neu- 
trally buoyant colloid drifting uniformly under an electric 
field in an unbounded fluid; and (II) a colloid with den- 
sity higher than the fluid, stably levitated by the balance 
between gravity sedimenting it towards a wall and an 



electric field driving it away. In I the excess noise vari- 
ance is 1 — 10 3 times the thermal noise variance; in II it 
is found to be at least as large as the thermal noise, and 
multiplicative in nature [see Fig.l (a)]. In both cases the 
extra noise is anisotropic - nearly an order of magnitude 
stronger along E than transverse to it. If S u and \u 
are the correlation and response functions of the colloid 
position at frequency u>, we see that u>S u ,/2lmXuj = T U1 
which should reduce to the thermodynamic temperature 
at thermal equilibrium, instead changes by a factor of 
2 as io changes from to 5.5Dk 2 (D = typical diffusiv- 
ity of the counterions and impurity ions hereafter collec- 
tively called microions, k _1 = Debye screening length), 
Fig 1(b). 




FIG. 1: (a) Plots showing the ratio of excess noise to thermal 
noise strength as a function of the distance from the wall z 
for different screening lengths. 

(b) Plots showing T w vs. ui for I and II where ui has been 
scaled by Dk 2 . 



We also find the steady-state solution P 00 (R) to the 
Fokker-Planck (FP) equation for the probability density 
of the colloid position R. We find that the effective po- 
tential U e ff(R) = -kBT\nPoo{R), where T is the ther- 
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modynamic temperature, is shallower than Squires's ef- 
fective potential for the same problem, Fig. 2. The ex- 
cess noise causes larger excursions than would arise from 
the bare thermal noise. Thus departures from equilib- 
rium behaviour should be evident even in single-particle 
experiments. We now show how these results were ob- 
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FIG. 2: The effective potential as computed from the Fokker 
Planck equation including the excess noise (solid line), com- 
pared to that in [6j (dotted line) for na — 4, a = 1/im, 0o = 4 
V for (a) iron core (red), (b) silica core. 



tained, beginning with case I. The colloid consists of 
a charge Q, surrounded by a counterion cloud of total 
charge — Q. The counterions plus ionic impurities, for 
short the microions, are assumed for simplicity to have 
identical mobility p and unit valency. For simplicity, we 
ignore advection of the microions by the velocity field. 
Thus the continuity equations for the densities n + and 
n~~ of positive and negative microions read 



at 



dn 
~dt 



= DV 2 n + - /LtV • (n + E) + V • (V2n+D/+), (1) 



D\7 2 n- + p\7 ■ (n~E) + V • (V2n~Df-), (2) 



where and / — are unit strength, independent white 
Gaussian noises, D = ksTp, is the ionic diffusivity, p 
being the Ohmic mobility, and E is the total electric 
field. 

Momentum conservation at zero Reynolds number, as 
appropriate for colloids, with thermal noise, leads to the 
fluctuating Stokes equation for the hydrodynamic veloc- 
ity field v of the incompressible suspension: 



where £ is a unit strength, divergence free conserving 
Gaussian white noise. We close the system with the Pois- 
son equation 

V-£=A (4) 

where the total charge density p — p® + p ml , p® = 
Q5{r — R(t)) and p m% = e(n + — ri~) being the colloid 
and microion charge densities respectively. The colloid 
is approximated as a point charge Q, and the effects due 
to nonzero size of the particle are taken care of by intro- 
ducing an ultraviolet cutoff 27r/a in Fourier space. We 
assume the colloid position R(t) simply moves with the 
fluid, 



8R(t) 
dt 



= v(R(t)). 



(5) 



The boundary conditions to be satisfied are: E{r — > 
oo) = Eq£ where Eq is the imposed field, and n ± (r — > 
oo) = no j the mean concentration of either species of ions. 
To make the problem tractable we linearize the equa- 
tions. The electrostatic force density in the Stokes equa- 
tion thus becomes p°EqZ where p° is the charge density in 
the absence of the electric field. The charge density has a 
steady and a fluctuating part: p°(r, t) — p(r) + 5p(r, t). 
Solving (|1I2[) for the charge densities in the absence of 
a field we get Sp(r,t) oc -i(q.f+, - q-f q J/(q 2 - i% + 
(ko) 2 ). The term 5pE z becomes an additive nonequi- 
librium Gaussian noise in the Stokes equation and thus 
contributes a noise proportional to Eq in the equation 
of motion for the colloid. Inserting the solution into the 
Stokes equation and evaluating v(R(t)), we find the ef- 
fective Langevin equation 



dR . . 



(6) 



for the colloid, where / is the thermal Gaussian white 
noise with variance 2fcsT/r, r cx i]a is the Stokes drag 
coefficient, C, with variance cx E 2 is the nonequilibrium 
excess noise and v — aQEon^ 1 /rja 2 is the steady elec- 
trophoretic velocity, a being a numerical factor of order 
unity. We measure the relative zero frequency strengths 
of excess and thermal noise by 



/ °°dt(Ci(o)0(t)) 
§/ °°di<A(o)A(t)) 



(7) 



and find a strong anisotropy in the excess noise and hence 
in 1Z; lZ Z z/7lxx — 8. Further, to estimate the typical 
scale of IZij, take (na) 2 = 2a 2 e 2 no/ekBT, and the ionic 



mobility p ~ e/rja m i c where 
then find that 



is the ionic radius. We 



77V 2 ?; - Vp + pE + ^j2k B Tr]i = 0; V 



0, 



(3) 



Tlzz oc 



eE 2 K 1 aa r , 



k B T 



(8) 
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which is the ratio of electrostatic to thermal energies 
in a volume involving three different lengths. Choos- 
ing typical experimental parameters E = 10 5 V/m, 
D = 10 _9 m 2 s _1 , na = 0.3, 1Z ZZ varies from 300 to 0.5 
for a varying from 10/im to 0.1/im. 

The excess noise should not be thought of as an 
effective temperature: it is not only anisotropic but 
also strongly frequency dependent. This is clearly seen 
through the ratio T u = ujS uj /21irXu J where S u and %w 
are the correlation and response functions of the colloid 
position as a function of frequency u>, evaluated as fol- 
lows. A body force density [-KR{t) + 5h(t)]6(r - R(t)) 
is added to the Stokes equation. The term in K holds 
the co-ordinate R at a stationary mean value and 8h(t) 
is a perturbation which leads to a shift 5R(t) in the 
particle position. We calculate \^ = (SR W )/Sh u and 
S U = J dte^iSR^SRit)) and find 



T 



1 + A 



1 + A 



1 



((q 2 + M 2 ) 2 + 



-d A q 



D 2 



V2(Kfl) 3 



M 4 + gr + 



leading to the results in Fig. lb. 

A more experimentally accessible system is case II, a 
colloid whose gravitational sedimentation towards a wall 
is countered by an electric field pushing it away. For 
this system the boundary conditions are <j)(z — 0) = 
</>o, <f>(z — > oo) = 0, and no slip and no penetration: 
v(x,y,0) = 0. 

The electrostatic boundary conditions imply that the 
mean microion density p av and electric field E vary with 
z. Since the fluctuating Maxwell stresses originate from 
the force density pE, the excess noise is now multiplica- 
tive, with a z dependence. The microions screen the po- 
tential so that the electric field decays away from the 
wall. Upto linear order in 4>o, ignoring the colloid charge 
and the microion density fluctuations, the screened elec- 
tric field is given by the linearized Poisson-Boltzmann 
equation solved with the above boundary conditions; 
E s = zE s0 e~ KZ , where E sn — fa/n -1 . For small col- 
loid charge Q and to lowest order in charge density fluc- 
tuations the force density in the Stokes equation is ap- 
proximately p°E s <d(z), where is the Heaviside func- 
tion and p° is the microion distribution in case I. The 
screened electric field recieves corrections due to p° and 
image charges required to satisfy boundary conditions, 
but for large screening and small Q this effect is seen to 
be very small. The buoyant weight of the particle adds 
a force density WS(r - R(t)) to ©. 

We use the Green's function Q(r,r') for the Stokes 
equation with a no-slip no-penetration wall in the xy 
plane as obtained by Blake [irj to solve ^ for the veloc- 
ity field v, including contributions from equilibrium ther- 
mal fluctuations and fluctuating Maxwell stresses pro- 



portional to E. Evaluating v at the colloid position R(t) 
gives the Langevin equation for the colloid in the form 
([5]), with an R dependent velocity Vq(R) with a zero 
corresponding to the minimum of Squires's effective po- 
tential. The thermal noise F and excess noise £ are now 
multiplicative, with correlations 

(C i (t,R(t))C j (t',R(t'))) = |d 3 fc^l F e-^ 2 + fe2 )( t - i ')w i (fe,^) ) 

x Hj(-k±,z(t'),-k z ), 
(F!(t,R(t))F;(t',R(t'))) = 5(t-t') [d 3 kn^(kx,z(t),k) 

x H' J (-k ± ,z(t),-k)dp, ( 

where the function Hi and H! i are the Fourier transforms 
of Qzi(r±, z, z') exp(— kz') and Gzi(r±, z, z') respectively 
with respect to r± and z' . By analogy to (7), we estimate 
the ratio of excess noise strength to the thermal noise, 
now as a function of z, by replacing z(t) and z(t') in ([9]) 
by a fixed z. 

The boundary conditions in the velocity field imply 
1Z F (0) = = 1Z^(Q); 1Z F rises to its saturation value 
within 1 — 2 particle radii, whereas TiS peaks at z = kT 1 
and falls off due to exponential fall in the field strength 
with increase in z so that the tail penetrates more into 
the bulk for smaller na. The peaks get sharper and the 
peak value larger with increase in na. Since n< and TZ F 
are anisotropic to varying degree, IZij is anisotropic and 
1Z zz /1Zxx is z dependent close to the wall saturating to a 
constant value close to 8, a few particle radii away from 
2 = 0. IZij scales as eVQKaa m i C /kBT; typically in ex- 
periments, Vq ~ a few volts, a = O.lpm to 10pm and 
Ka r~j 1 so that the ratio is at least 0(1) within a few a of 
z = k^ 1 . The ratio has dipped in magnitude compared 
to case I due to exponential fall in the field magnitude 
and the wall boundary conditions. 

As a result of the excess noise proportional to E, the 
steady-state probability distribution P x (R) is not given 
by a Boltzmann weight determined by Squires's effec- 
tive potential. To obtain we must solve the equation 
for R treating carefully the multiplicative and coloured 
noises in (6). The equation is driven by a continu- 
ous, weighted, linear superposition of independent, mul- 
tiplicative Ornstein-Uhlenbeck (OU) [13j noise sources. 
We construct the Fokker Planck equation corresponding 
to (6) by generalizing the results of [3] for a single mul- 
tiplicative OU noise. Since the weights and relaxation 
times of the constituent noises in our problem are finite 
as a result of Debye screening, all integrals arising in this 
procedure are convergent. The resulting FP equation 
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takes the form 

d{WP) 



dP 
~8t 



dz 



oz oz 



dT—[G T — (G T + tG t W - tG' t W)P] 
oz oz 



(12) 



where we have schematically replaced integrals over wave 
number by integrals over relaxation time r; r goes as 
[D(k 2 + k 2 )}- 1 and G T {R) = (k 2 + K 2 ) _1 W(p, z, k). 
Poo(R) is obtained by solving (TH?)) for the zero flux 
condition and U e ft = —kBT\nP 00 (R). The probabil- 
ity distribution P^ analogous to that of contain- 
ing only the bare thermal noise, is obtained by setting 
G T = in (fl~2f : the corresponding effective potential 
E/Pjj = — ksT In P^. The two potentials show large dif- 
ferences near the wall where the excess noise is signifi- 
cant. S R = {U° ff )" \ Rav /(U e ff)" \ Rav , gives the relative 
confining strengths of the two potentials and Sr » 1 
for R av close to the wall (Sr = 2.7 for Fig. 2.); i.e. U e ff 
is shallower than U^j. The excess noise causes the col- 
loid to explore a wider range of z. This nonequilibrium 
effect is strongest if the colloid is heavy, so that its mean 
position is close to the wall. Fig. 2 shows the comparison 
between U e ff and U®ff for a colloid with mass density 
corresponding to (a) silica and (b) iron. 

Lastly, let us check that the neglect of the advection of 
microions by the hydrodynamic velocity field v was not 
too bad. The appropriate Peclet number Pe= v/Dn 2 a, 
the ratio of the rate at which the colloid shears the 
medium to the rate at which microion densities relax. 
The Langevin equation for the colloid implies a typical 
deviation from mean position and hence a typical value 
for the right-hand-side of (6). Using this, with a sur- 



face charge density 



10' 



10 3 C m 2 which is 



reasonable approximation for silica particles 11| , we find 
Pe= 0.35, 0.12, for na — 3, 5 respectively. Thus the 
neglect of advection is not a bad approximation. 

In summary, we have formulated the statistical dy- 
namics of a single colloidal particle in a static electric 
field. We show that thermal agitation of counterion and 
impurity charge densities, in the presence of the im- 
posed field, leads to noisy Maxwell stresses and hence 
to an additional noise term, proportional to the field, 
in the effective Langevin equation for the colloidal par- 



ticle. This noise is nonequilibrium in nature and highly 
anisotropic, and leads to strong, frequency-dependent de- 
partures from the Fluctuation-Dissipation Theorem. The 
effective potential as inferred Q by combining Stokes 
drag with the mean velocity when the particle is dis- 
placed from its steady-state position differs substantially 
from that obtained by taking the logarithm of the steady- 
state probability distribution. Our results are quantita- 
tively testable, e.g., in experiments such as those of 0, 
and should form an essential ingredient in understanding 
the nonequilibrium steady states of colloids in electric 
fields. Generalizations of our treatment to include oscil- 
latory fields, as well as more than one colloidal particle, 
are underway. 
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